Nonequilibrium probabilistic dynamics of the logistic map at the edge of chaos 
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We consider nonequilibrium probabilistic dynamics in logistic-like maps x t +i = 1 — a|a;t| z , (z > 1) 
at their chaos threshold: We first introduce many initial conditions within one among W » 
1 intervals partitioning the phase space and focus on the unique value q sen < 1 for which the 



entropic form S q = 
vanishes like t~ 1 ^ qrc 



(W)- 



increases with time. We then verify that S qsen (t) — S qsen (oo) 



[q r ei(W) > 1]. We finally exhibit a new finite-size scaling, q re i(oo) — 
qreiiW) oc VK~^ ae "'. This establishes quantitatively, for the first time, a long pursued relation 
between sensitivity to the initial conditions and relaxation, concepts which play central roles in 
nonextensive statistical mechanics. 
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Connections between dynamics and thermodynamics 
are not at all completely elucidated. Frequently ther- 
mostatistics may sound as if it would consist only of a 
self-referred body, which could dispense dynamics from 
its formulation. But this has long been known not to be 
true (see [1] and references therein). One possible rea- 
son for this essential point having been poorly empha- 
sized is that when dealing with weakly interacting sys- 
tems, Boltzmann-Gibbs (BG) thermodynamical equilib- 
rium may be formulated without referring to the under- 
lying dynamics of its constituents. As complex systems 
came to the front line of research (by complex we mean 
here the presence of at least one of the features: long- 
range interparticle interactions, long-term memory, frac- 
tal nature of a pertinent phase-space, small-world net- 
working) , it became necessary to revisit this fundamental 
issue [2]. Indeed, a significant number of systems, e.g., 
turbulent fluids [3,4], electron-positron annihilation [5], 
economics [6], motion of Hydra viridissima [7], and oth- 
ers, are known nowadays to be not properly described 
with BG statistical mechanical concepts. Systems such 
as these have been successfully handled within nonexten- 
sive statistical mechanics. The purpose of this Letter is 
to numerically illustrate, with a very simple dynamical 
model (logisticlike maps), that basic dynamical concepts 
such as the sensitivity to the initial conditions and the re- 
laxation towards equilibrium are deeply entangled, thus 
yielding a new analytic scaling connection. 

The basic building block of nonextensive statistical me- 
chanics is the nonextensive entropy [8] 
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2-~n=l Pi 



( q e 11). 



(1) 



The entropic index q characterizes the statistics we are 



dealing with; q — 1 recovers the usual BG expression, 
Si = — Y^iLi Pi m Pi- We may think of q as a biasing 
parameter: q < 1 privileges rare events, while q > 1 
privileges ordinary events: p < 1 raised to a power q < 1 
yields a value greater than p, and the relative increase 
P q Ip = P q ~ 1 is a decreasing function of p, i.e., values of p 
closer to (rare events) are benefited. Correspondingly, 
for q > 1, values of p closer to 1 (ordinary events) are 
privileged. BG (q = 1) is the unbiased statistics. A con- 
crete consequence of this is that the BG formalism yields 
exponential equilibrium distributions, whereas nonexten- 
sive statistics yields power-law distributions (the BG ex- 
ponential is recovered as a limiting case: we are talking of 
a generalization, not an alternative) . Many developments 
concerning this formalism are available in the literature, 
[9,10], but there are still important points that remain 
to be better understood. One of them is what lays be- 
hind the quite impressive results recently put forward for 
fully developed turbulence. Agreement between theoret- 
ical treatments and experimental data has been achieved 
along two different approaches, both based on nonexten- 
sive statistical mechanics, but, at first sight surprisingly, 
one of them with q > 1 [3] and the other with q < 1 
[4]. A similar scenario (using both q < 1 and q > 1) 
is also found in the description of low dimensional dy- 
namical systems, particularly in the z-logistic maps at 
the threshold of chaos. This Letter intends to shed light 
on this problem, by investigating these maps both in the 
chaotic region and at the edge of chaos, and by showing, 
for the first time, a long pursued connection, namely, 
between the sensitivity-based and the relaxation-based 
approaches. 

Consider the z- logistic iterative equation Xt+\ = 1 — 
a\x t \ z (-1 < x t < 1; < a < 2; z > 1; t = 0,1,2...) 
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for which z = 2 recovers the logistic map, and parti- 
tion the phase space into W cells of equal measure, pi 
being the probability to be soon associated with the i th 
cell. Several works [11-15] have shown that the descrip- 
tion of such dynamical systems is properly achieved by 
using the nonextensive entropy, Eq. (1), with a specific 
^-dependent value of q. 

Let us first recall the sensitivity-based approaches, 
which determine q sen < 1 (sen stands for sensitivity). 
Up to now, three different methods have been devel- 
oped which provide the special value q sen (sometimes 
noted as q*). The first method [11] is based on the 
sensitivity to initial conditions. At the threshold of 
chaos [a — a c (z)}, the Lyapunov exponent Ai vanishes 
and the system is weakly sensitive (weak chaos), i.e., 
the upper bound separation between nearby trajectories 
£(t) = lim Aa; (o)— o 2^ typically evolves [16] in time ac- 
cording to the power-law 

£(*) = [l + (l-«aen)A,. aB t]l=£^ (A 9se „ > 0), (2) 

which is solution of the nonlinear differential equation 
£ = \ scn £, qs "" (this equation recovers, for q sen = 1, 
the usual exponential sensitivity to the initial condi- 
tions, referred to as strong chaos). The second method 
is based on the geometrical description of the multifrac- 
tal attractor [12,13]. The value of q sen is determined by 

t— ; — = 7T— i where a min and a max are the val- 

ues of the end points of the multifractal function /(a). 
This beautiful equation relates dynamics (left-hand side) 
with geometry (right-hand side). The third method is 
related, as we shall specify, to a conjectured general- 
ization of the Pcsin relation [14,15]. In one of the W 
cells in which the phase space has been divided, we ini- 
tially put (randomly or uniformly distributed) N points 
(N ^> W). We then follow the occupancies {Ni(t)} of all 
cells Ei^li-^W — which enable the calculation of 
the probabilities Pi(t) = Ni(t)/N, which enable in turn 
the calculation of S q (t) as given by Eq. (1). We can next 
define a nonextensive version of the Kolmogorov-Sinai 
entropy, namely K q = lim^^ limvy^oc limAr^oo S q (t)/t. 
The special value q — q sen is the one for which the en- 
tropy production K q is finite (if q < q sen , K q — > oo, and 
if q > q sen , K q — ► 0). It is quite remarkable that all three 
methods give (within an acceptable numerical error) one 
and the same value for q sen (q sen = 0.2445... for z = 2). 

Let us now turn to the relaxation-based approach of 
the problem, first tackled in [17]. We start now [17] with 
an ensemble of initial conditions uniformly spread over 
the entire phase space (a procedure that resembles the 
Gibbsian approach of statistical mechanics), instead of 
within a single among the W cells (closer to a Boltzman- 
nian approach), and investigate the rate of convergence 
onto the multifractal attractor at the onset of chaos. At 
t = all W cells are occupied (hence S q — w 
As t increases, the number W occup i e d of cells that have at 



least one point inside typically decreases (shrinking of the 
Lebesgue measure) as 1/[1 + (q re i — 1) t/T q ] 1 /^ qrel ~ v> (rel 
stands for relaxation; q re i > 1; r q > is a characteristic 
relaxation time). 

When applied to the chaos region (e.g., z = 2, a = 2), 
all these methods yield q sen = q re i = 1, in accor- 
dance with the usual BG concepts, founded on strongly 
chaotic systems, consistently with Boltzmann's "molec- 
ular chaos" hypothesis (see [18] and references therein). 
When applied to the edge of chaos a c (z) of the z-logistic 
maps, it is obtained [17] q re i(oo) [g re ;(oo) = 2.4... for 
z = 2]. In this Letter, we start with N points (we adopt 
N = 10W for numerical convenience) uniformly dis- 
tributed inside a single specific cell chosen as soon spec- 
ified. We compute S qsen (t) and verify that, for a fixed 
value of W, S qsen (t) asymptotically reaches S qsen (oo) 
[S qeen (oo) > monotonically increases with W\. S qsen (t) 
reaches its t — > oo value from above, preceded by an over- 
shooting, which might be very strong. The initial cell is 
chosen as that one which presents the highest overshoot- 
ing. 

If we consider the strongly chaotic case a = 2, hence 
q sen — 1, S\(t) slightly overshoots and rapidly saturates, 
as illustrated in Fig. 1(a) for z = 2. At the onset of chaos, 
however, S qgen (t) presents a very pronounced overshoot, 
and slowly approaches its final value S qsen (oo) [Fig. 1(b)]. 
We then compute, for all times after the overshooting, 
AS qsrn (t) = [S q3en (t) - ^^(oo)]. We find a power-law 
decay, whose log-log slope depends on W. (Log-periodic 
oscillations such as those reported in [17] are also de- 
tected here, and are better visualized for z < 2.). We 
identify this slope with \/(q re i(W) — 1) (as an alterna- 
tive way of obtaining q re i). In the limit W — > oo, we 
observe that q re i(oo) precisely coincides, for all values of 
z, with the value found in Ref. [17] (their values corre- 
spond to infinitely fine graining W — > oo). The power- 
law region of AS qsen (t) increases with W (see Fig. 2). 
In order to better identify this region, we consider only 
those time intervals whose linear correlation coefficient 
is larger than 0.99. In order to increase the precision of 
qrei (W) we take averages over various possible intervals, 
more precisely, ranging from various initial and various 
final times, as shown in Fig. 2. Moreover, in order to fur- 
ther minimize the effect [on the precision of q re i(W)\ of 
the oscillations of S qeen (t) at the onset of chaos, and also 
to better identify the attractor — the invariant distribu- 
tion — we introduce a numerical improvement in the iter- 
ation rule which does not affect the asymptotic dynamics 
of the map: we average the distribution corresponding to 
x t with that corresponding to x t -i (the actual values of 
x must be preserved while averaging). This procedure 
strongly reduces the considerable fluctuations in the dis- 
tribution of points at t — > oo and enables us to identify 
the attractor associated with the map (we refer to the 
attractor in the space of distributions, and not in phase 
space). The computer implementation of these numeri- 
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cal benefits demands additional memory and an increase 
of the CPU time. The oscillations that remain are con- 
siderably smaller, and they can be further reduced by 
extending the method in order to take into account more 
than one previous time step. However, for the accuracy 
we seek in this work, it is clearly enough to average the 
distributions of x t -\ and xt- We emphasize that this kind 
of averaging is just a trick to minimize the fluctuations 
(the original dynamics is preserved, once we average the 
distribution of x, and not x itself). If the procedure is 
not adopted, we find the same results, though with less 
precision. After all these numerical improvements have 
been implemented, a remarkable law emerges, namely, 

Qrel(oo) - Qrel(W) OC tU~' 9sen ' (3) 

(see Fig. 3). This equation is somehow analogous to the 
expressions that appear in finite-size scaling [19]. This 
law is quite impressive, since it relates two basic dynam- 
ical properties of dissipative systems, namely, relaxation 
onto equilibrium and sensitivity to the initial conditions 
(the basis for mixing). Moreover, the (coarser or finer) 
graining - l/W for the z-logistic maps — is involved, 
pretty much in the same way as it occurs with the values 
of q (here noted as q re i) appearing in Beck's approach 
[3] of fully developed turbulence, where the entropic in- 
dex depends on the distance between the two points at 
which the fluid velocities are measured. Other phenom- 
ena where a similar scaling relation may occur are the 
distribution of transverse momenta of hadronic jets pro- 
duced in electron-positron annihilation experiments [5], 
and the saddle-point dynamics of the Henon-Heiles sys- 
tem [20]. Let us also mention that the scaling relation, 
Eq. (3), has a moderate sensitivity to the value used 
for q sen (see the Inset of Fig. 3). It is important to 
emphasize that Eq. (3) reflects the fact that the lim- 
its lim t ^oo lim^/^oo (relevant for q sen ; see for instance 
[14]) and lim^^co lim^oo (relevant for q re i, as exhib- 
ited here) do not commute (nonuniform convergence) in 
general (in other words, generically q sen ^ q re i, whereas 
for strong chaos it is q sen — q re i = 1). This result 
is similar to the noncommuting t — > oo and M — > oo 
limits for long-range interacting M-body Hamiltonian 
systems [21], for which vanishing Lyapunov exponents 
have already been detected [22] (we remind that these 
limits do commute for short-range interactions). We fi- 
nally exhibit Fig. 4, which shows q re i(oo), as well as q sen 
(taken from [11]), for typical values of z. For q re i(oo) 
we have indeed verified the coincidence (with high preci- 
sion for all z) with the values observed in [17] by starting 
with uniform occupation of the phase space and follow- 
ing the shrinking of the Lebesgue measure. In addition 
to these results, we observe that q sen {z — > oo) ~ 0.72 and 
limw^oc q re i(z -> 1) ~ 4/3. 

In conclusion, the scenario which emerges is that the 
thermal equilibrium of all Hamiltonian systems whose 



dynamics is dictated by noninteracting or short-range 
interacting particles is described by Boltzmann-Gibbs 
statistical mechanics, and usual thermodynamics apply, 
typically together with the usual exponential relaxation 
to equilibrium. But for systems with complex dynamics 
(complex in the sense previously described), we cannot 
know a priori how the relaxation to the stationary state 
occurs. However, for a vast class of such systems, power- 
law relaxation does occur. To determine the associated 
exponent we must analyze at least once the dynamics of 
that particular nonextensivity universality class, in such 
a way as to determine once the corresponding values of q. 
For z-logistic maps, Eq. (3) (main result of the present 
work) gives the connection between those values of q. 
The analogous task for Hamiltonian systems would be 
very welcome. 
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FIG. 1. S qssn (t) for z = 2 . (a) Chaotic region, with q sen = 1; (b) Edge of chaos, with q sen — 0.2445, associated to the 
region before the peak [12,14], and with q re i{W) > 1, associated to the region after the peak, (the highest value attained by 
So. 2445, for a given W , is about 70% of the value corresponding to equiprobability) . Notice that for a = 2, the stationary state 
is achieved for t much smaller than that for a — a c , which of course reflects the fact that relaxation is exponentially quick in 
the former (as it is also the case for a < a c ) whereas it is a power-law in the latter. 
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FIG. 2. AS qaen (t) = Sq acn (t) — S qaen (oo) at the 2 = 2 chaos threshold, for typical W. For each W we consider as the 
power-law region (linear correlation coefficient above 0.99) that going from the left interval to the right one. Furthermore, the 
values for q re l{W) in Fig. 3 were obtained by averaging over many starting (ending) points inside the left (right) interval. 
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FIG. 3. IV-dependence of q re i for typical values of z. The W — > oo extrapolated values coincide with the values reported in 
[17]. The abscissa has been chosen to be (W/IOOO) - ' 1 '"" (instead of W~' 9se ™') for better visualization. Inset: Linear correlation 
coefficient procedure which has been used to select a specific numerical value for the sensitivity entropic index q ae n{z). 
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FIG. 4. z-dependences of q se n (from [11], or self-consistently through Eq. (3)) and q re i(oo). Inset (a): z — » 1 extrapo- 
lation for g re ;(oo); q re i(oo) — 4/3 + 1.077(2 — l) 058 fits the z < 3.0 data. Inset (b): z — > oo extrapolation for q sen (oo); 
q S en{z) = 0.72 - 1.525/ z 1 ' 7 fits the z > 1.75 data. 
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